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Results from an extensive relativistic many-body analysis utilizing a realistic effective QCD Hamil- 
tonian are presented for the meson spectrum. A comparative numerical study of the BCS, TDA 
and RPA treatments provides new, significant insight into the condensate structure of the vacuum, 
the chiral symmetry governance of the pion and the meson spin, orbital and flavor mass splitting 
contributions. In contrast to a previous glueball application, substantial quantitative differences are 
computed between TDA and RPA for the light quark sector with the pion emerging as a Goldstone 
boson only in the RPA. 
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T-H ^ Common to the diverse areas of condensed matter, molecular, atomic and nuclear physics is the routine implemen- 

^ . tation of many-body techniques such as the Bardeen, Cooper, SchriefFer (BCS), Tamm-DancofF (TDA) and Random 

^ Phase Approximation (RPA) methods. Particle physics, with an inherent few-body nature, has generally been devoid 

^~~i of such applications even though hadronic structure, requiring a relativistic QCD description, is an extremely challeng- 

^^ ing many-body problem. The purpose of the present letter is to report a comparative study documenting the powerful 

^~^ . utility of the above techniques for hadronic systems and to detail new, important meson structure results clarifying the 

nature of spin splittings and role of chiral symmetry. The equations of motion, while numerically solvable, exhibit a 

K^ . richness and complexity beyond the simple two-body equations such as the Bcthe-Salpeter or generalized Schrodinger 

^s. I schemes. We find that both TDA and RPA solutions to an approximate QCD Hamiltonian with linear confinement 

l/^ . reproduce the meson spectrum except for the pion, where only the RPA reasonably describes the mass and decay 

C^ ' constant due to proper implementation of chiral symmetry. 

This work complements our previous many-body treatment pi of the giuonic sector in which the lattice gauge 
"measurements" were reproduced. Our collaborative program seeks to develop a rigorous effective Hamiltonian from 
Qs i QCD and then to comprehensively investigate hadronic structure by systematic, accurate diagonalization utilizing 
controllable approximations. Ref. [^ details our renormalization program, based upon a continuous cut-off regular- 
ization and similarity transformation. That work addressed only the gluon sector but similar effort is currently in 
1 '' progress for the quark sector. Accordingly, this paper presents many-body solutions for only the unrenormalized 
^H, effective Hamiltonian. The starting point is the approximate QCD quark Hamiltonian in the Coulomb gauge 

-: H = j dx^\{x){^ia-V + Pm,)^>,{x)-]^j dxdyp%x)V{\x^y\)p-{y), (1) 

/^ ' involving the quark field ^q{x), current quark mass ruq and color density p°'{x) = '^Ux)T°''^q{x). Coupling to the 
giuonic sector is omitted and the Faddeev-Popov determinant is replaced by its lowest order unit value. Consistent 
with our previous work [|l|, the confining potential is a linear interaction, V = a\x — y\, rather than the harmonic 
oscillator g,^ since lattice gauge theory generates this form with slope (string tension) a — 0.18 GeV^. Instead 
of a simplified gap differential equation for the harmonic oscillator potential, we solve a numerically quite sensitive 
nonlinear integral equation (see below) and reproduce earlier results |g]. The density-density two-body form permits 
only color singlets in the physical spectrum as other SUc{S) representations are shifted to infinite energy. 

Next we introduce our first many-body improvement by performing a BCS rotation (similarity transformation) 
from the bare (undressed) quark basis to an improved quasi-particlc basis. This entails rotated spinors in the quark 
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field expansion with quasi-particle operators B, D instead of bare operators b, d (see ref. M for details) 

*«=e/(|^ h(^)^^(^)+^>(-^)^i(-^)]^^'''^- (3) 

The spin state is denoted by A and color indices are suppressed. The gap angle (j){k) governs the BCS vacuum, |rj), 
defined by Bx\^) — D\\^) = 0. This vacuum, a coherent state containing quark condensates (Cooper pairs), is an 



improvement over the trivial vacuum. The gap angle is obtained variationally by minimizing the vacuum ground 
state 6{fl\H — E\^) = yielding the gap equation which is similar to the Schwinger-Dyson equation for the quark 
self-energy in the rainbow approximation 



2 f dp 



ksk-mqCk = - / , ,^ V{\k~p\)[skCpk -p- SpCk], (4) 



wheieV{\k—p\) — ~8Tra/\k~p\ is the linear potential in momentum space. The solution Sfc = sin<^(A:), Ck = cos0(fc) 
also provides the quark condensate (qq) 

(qq) = {n\^iO)^{0)\n) = -^ fdpp^ Sp, (5) 



by Sfc = TOj?*'"(A:)/v/mg^"(fc)2 + fc^ and obtain rhq from the slope of the gap angle near zero momentum (m, 



in the BCS vacuum. From the gap equation at large k, Sk — > nig / y/rriq^+k^ , yielding a quadratically divergent 
condensate for non-zero current quark mass which must be renormalized 0. For nig = we compute {qq) ~ 
— (113 AleV)^. We also added the Coulomb interaction with a reasonable cutoff and found a slight improvement to 
119, in agreement with the more elaborate, renormalized result of ra] but still substantially less than lattice theory 
results (« 250). The BCS vacuum also exhibits spontaneous chiral symmetry breaking resulting in a constituent 
quark mass, rhq, which can not be extracted from the gap energy, e^, 

ek ^ niqSk + kck ~ - / .„ ^-, ^{1^ ~ p\)[cpCkk ■ p + SpSk], (6) 
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since Cfe diverges for fc — > 0. Following previous prescriptions &-H], we introduce the running dynamic mass, TOj?*'"(fc), 

q ~ 

m^''"(0)). This yields m^/d ~ 80 MeV for u, d flavors with m„ — rrid — 5 MeV and rhs « 250 MeV for the s quark 
with TTis = 150 MeV. Again these values are somewhat lower than used in phenomenological quark models indicating 
a more sophisticated vacuum is needed which is provided by the RPA as detailed below. 

Using these quasi-particle creation operators we now address excited meson states and first construct the TDA 
Fock space wavefunction built on the BCS vacuum. For a meson with quantum numbers nJ'^ (radial/node number 
n, total angular momentum J and parity tt) the leading qq state is given by 

LSmLtns 

Diagonalizing H in this model space generates the TDA equation of motion (analogous to the Bethe-Salpeter equation) 

{^-'^\[H,BiDl]\n) = {E„j. E,)^-';. (9) 

Very significantly, the relativistic structure of our effective interaction contains an important spin dependence. This is 
revealed more clearly in the TDA partial-wave equations for a meson in state nJ'" with mass Mnj^ having quantum 
numbers L (orbital) and S (total spin) 
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(M„j, - 2ek)rLi ik)^J ^^ Kis(k,p) rd (P), (10) 



with kernel ii'/g(fc,p) structure: 1) pseudoscalar (J'^ = ), 2{ckCpVi + {1 + SkSp)Vo); 2) scalar (J^ = 0+), 2{ckCpVo + 

(1 + SkSp)Vi); 3) vector (J'^ = 1^), 2ckCpVi + {1 + Sfc)(l + Sp)Vo + (1 - Sp){l - Sfe)(|f2 - |K)); 4) pseudovector 

{J^ = 1^), CfcCp(Vo + V2) + 2(1 + SpSk)Vi. Here Vi is the angular integral over x — k ■ q oi V{\k — q\) with powers a;\ 

The TDA spectrum is given in Fig. 1 (dotted lines) for the pseudoscalar and vector mesons. Note that in this model 

the p and uj states are degenerate and there is no dynamic flavor mixing (the standard SUf(S) flavor mixings have 

been adopted for rj, 77', to and 4>)- Considering that a is the only, but fixed, parameter in this model, the spectrum 



is qualitatively reasonable with the exception of the ground state pion. From model calculation the spin splitting 
between the p and tt is about 200 MeV , insufficient to describe the roughly 600 MeV observed difference 0. This 
shortcoming is due to the inability of the TDA to properly include constraints from chiral symmetry. 

Finally, we formulate the RPA |^ and in analogy to nuclear physics applications generalize the meson creation 
operator, Q\j„ , for state nJ'^ to both create and destroy qq pairs with an improved vacuum satisfying Qnj^ \^rpa) = 
and containing quark correlations beyond the BCS. To obtain the RPA equations of motion we make use of the quasi- 
boson approximation in which pairs of operators, BD, are treated as boson operators. For the important pseudoscalar 
(J'^ = 0~) meson channel we obtain the coupled partial- wave RPA equations 

1 r°° dpj^ 
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2efeX"(fc) + - / ^ [X^{p)F{k,p) + r"(p)G(fc,p)] = M„X"(fc), 



1 f°° dv T? 
2ekY^{k) + - / -P^ [Y^{p)F{k,p) + X^{p)G{k,p)] = -AUY^{k), 



3^0 (2^)' 
where 

F{k,p) = 2cpCkVi+2{l + SpSk)Vo, (11) 

G{k,p)^2cpCkVi-2{l~SpSk)Vo, (12) 

with similar expressions for the other meson channels. The RPA spectrum (dashed lines) is summarized in Fig. 
1. Notice the improvement with only the pseudoscalar states (tt and 77) shifted downward from the TDA. Related, 
we now also obtain the correct chiral limit for the pion mass. As expected from chiral arguments concerning the 
Goldstone boson nature of the pion, the RPA pseudoscalar mass (pure u/d flavor) approaches zero for niq ^ as 
illustrated in Fig. 2 (solid curve). Appropriately, the scalar meson, /o, mass (dotted curve) converges to a non-zero 
value in the same limit (the scalar and pseudovector meson spectrums are also reasonable but not shown due to space 
limitations). Hence the RPA clarifies the major source of mass splitting between the tt and p - roughly 200 MeV from 
spin dependence, as in the TDA, but a much larger amount, about 400 MeV , due to chiral symmetry constraints. 

The RPA also provides impressive improvement in the quark condensate and pion decay constant. Performing an 
expansion of the RPA vacuum in powers of boson operators and keeping only leading corrections exciting up to two 
mesons from the BCS vacuum, we obtain a much larger condensate. For zero current quark mass, (gq) ~ — (300Af e^)'^, 
in better agreement with the currently accepted value. Interestingly, this overestimation of ground state correlations 
seems to be characteristic of RPA, as previously documented in other fields of physics |g| . Finally, in the chiral limit 
we compute the pion decay constant /^ = 120 AleV , in reasonable agreement with data (93 MeV) and substantially 
better than the TDA value of 10 MeV . Related and also noteworthy, we have numerically verified the generalized 
Gell-Mann-Oakes-Renner relation 



2m,{qq) ^YNllfl (13) 



by independently computing terms on both sides of the equation. The excited pseudoscalar states negligibly contribute 
as their /„ are suppressed in the chiral limit. 

Summarizing, we have documented the utility of several many-body techniques, especially the RPA, for investigating 
the QCD structure of hadrons. In conjunction with our previous glueball study we have also further established our ef- 
fective Hamiltonian and many-body approach as, with just the independent lattice parameter a, the semi-quantitative 
features of the vacuum, meson and glueball spectrums have all been reproduced. Because the approximations are 
controllable, this framework is amenable to systematic improvement and should be appropriate for more challenging 
hadronic investigations, some in progress, such as baryons, hadron hidden flavor (higher Fock states) and hybrids. 
The details of the present calculation, as well as several issues concerning renormalization [^ and the rj ~ rj' system 
are deferred to a subsequent, major publication. The authors wish to recognize useful discussions with J. E. Ribeiro, 
A. P. Szczepaniak and the NCSU theory group. F. L-E. acknowledges a SURA- Jefferson Lab fellowship. This work 
is partially supported by grants DOE DE-FG02-97ER41048 and NSF INT-9807009. 
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FIG. 1. Pseudoscalar and vector meson spectrum. Data (solid), RPA (dashed) and TDA(dottcd). 
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FIG. 2. Chiral symmetry in the RPA. For m, -^ the pseudoscalar (sohd) but not scalar (dotted) meson mass vanishes. 



